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Abstract. The evaluation of iterated primitives of powers of logarithms is 
expressed in closed form. The expressions contain polynomials with coeffi- 
cients given in terms of the harmonic numbers and their generalizations. The 
logconcavity of these polynomials is established. 



1. Introduction 

The search for closed forms of definite integrals has the extra appeal of connecting 
many diverse areas of mathematics. At the beginning of Calculus, the student is 
usually told of 

/oo 
e~ x2 dx = v^F, 
-oo 

that hints to a relation between exponential and trigonometric functions. The 
reader will find in Chapter 8 of [3] a collection of different proofs of (II. ip . These 
include the classical ones, presented in most textbooks, as well as one by R. A. 
Kortram relating this evaluation to the number of representations of a number 
as a sum of squares. 

The second author has begun the project of establishing all the entries in the 
classical table of integrals by I. S. Gradshteyn and I. M. Ryzhik [§]. An example of 
the surprising connections encountered in this process is the subject of this note. 

Consider the sequence of functions {f n (x) : n € N} defined by the iterated 
integrals 

(1.2) f (x) = ln(l+x) 

fn(x) = f fn-l(t)dt. 
JO 

The first few examples are given by 

/i(x) = -x + (1 + x) ln(l + x) 

Mx) = -|(3x + 2) + i(l + .T) 2 ln(l + x) 

f 3 (x) = -|;(llx 2 + 15x + 6) + i(l + a; ) 3 ln(l + .T) 

h( x ) = -7^(25a; 3 + 52a; 2 +42.T + 12) + T i r (l + a;) 4 ln(l + a:). 
zoo 24 

This data suggests that 

(1.3) fn{x) = -xA n {x) + B n {x) ln(l + x) 
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where A n and B n are polynomials. Moreover, the value 



(1.4) B n {x) = -Al + x) n 

can be guessed from the data above. 

In view of the fact that a closed form for the polynomials A n {x) seems harder 
to find, we begin by considering 

(1.5) a n := lcm{ Denominators of A n (x)}, 

the sequence of denominators in the reduced form of A n (x). This sequence starts 
with 

(1.6) {1, 1, 4, 36, 288, 7200, 43200, 2116800, 33868800}. 

A direct search in Neil Sloane's [15] gives no information. On the other hand, the 
quotient 

(1-7) A 



n • 

na n -i 



gives the values 

(1.8) {1, 2, 3, 2, 5, 1, 7, 2, 3, 1, 11, 1, 13} 

and this is the sequence ^4014963 of Sloane. Namely 



(1.9) $ n 



p if n is a power of p 
1 otherwise. 



The sequence p n is the exponential of the von Mangoldt function 
(1.10) A„ = 

one of the basic functions in the theory of prime numbers [9]. 



lnp if n is a power of p 
otherwise, 



In this note we provide an explicit formula for the polynomial A n (x) in terms 
of harmonic numbers. This establishes the form of j3 n discussed above. Section [5] 
considers the iterated primitives of a power of ln(l + x). A sequence of polynomials 
involving the generalized harmonic numbers is given. 

2. The recurrences 
The first few values of f n (x) suggest the ansatz 

(2.1) fn(x) = -xA n (x) + B n (x) ln(l + x), 

with A n and B n polynomials in x. This is replaced in the relation f' n {x) = f n -i(x) 
to produce 

(2.2) B' n (x) = B n _ x {x) 

(2.3) xA' n (x) + A n (x) = x A n ^{x) + ^^. 

1 + x 

The expression for B n (x) in (|1.4[) is obtained directly from here. 



ITERATED PRIMITIVES OF LOGARITHMIC POWERS 



3 



Theorem 2.1. The polynomial A n is given for n > by 



(2.4) A n {x) = -^( n h ){H n -H n _ k )x 



fc=i 



fc-i 



where H n = l + ~ + j|+--- + iis the harmonic number. 



Proof. The formula is proved by induction. Write 



(2.5) A n (x)^J2 a ^ k ~ 1 

k=\ 



and use the recurrence (|2.3[) with the explicit expression for B n to conclude that 
a n,i = l/n\ and for 2 < k < n 



, . 1 Ai — 1 

(2.6) fca n , fe = a„_i ife _i + — 

rt! \ k — 1 



This recurrence now gives (|2.4j) . □ 



Note 2.2. Iterated integrals of ln(l — x) were considered by Mathar [T31 section 4.3], 
who obtained a similar recurrence. 

Theorem 2.3. The polynomials A n (x) are given by 



fc=0 v 7 m=l 
Proof. The definition of /„ (cc) implies that 

,2 - 8) /»M = ]C j(j + i)~ ■ (j + n) ' 

Using the partial fraction decomposition 

(29) I = L f(-ll(n 

j(j + !)••• (j + n) n\f^ o j + k\k 
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produces 

j=l k=0 J v 

= -*"gM)*(:)(-*>-'§^ 

fc=0 x 7 m=fc+l 

- -e(:W-m 1+ *>-x;^ 

fe=0 v 7 \ m=l 



(l + x)«ln(l + + 



fc=0 m=l 



k=o x 7 m=l 

Dividing this sum by — x gives the result. □ 

Comparing the two expressions produced for the polynomials A n yields the next 
identity. 

Corollary 2.4. Let neN. Then 

(2.10) 2J ( fc J 21 £ m "X/./J 21 X/jTj-j-n-fc" 

fc=0 v 7 m=l fe=l v 7 m=l 

3. The denominators of A n (x) 

In this section we analyze the polynomial A n (x) and compute its denominator 
when written in reduced form. The proof employs some elementary number theory. 

Theorem 3.1. For n > 1, the common denominator for A n (x) is given by 

(3.1) a n — n\ lem (1, 2, . . . , n). 

The proof of the theorem employs a preliminary divisibility result. 
Lemma 3.2. Let p be a prime, n > 1, and > such that p k < n. Then the 
power of p dividing the denominator of ^ (iJ„ — H n _ p k) is p fc . 

Proof. Observe that 

P fc -i x 

(3.2) H n -H n _ p k = \ :. 

i=0 

In (|3.2[) , the index i ranges over a set of length p k ; therefore there is a unique 
index i such that p k divides n — i , namely io = n — p k \n/p k \. For i ^ i write 
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n — i — p ai Pi and n — io — p a /3, with a > k, on < k — 1 and p not dividing j3, Pi 
Therefore 

1 1 

can be expressed as 

TT -H - + J_ 

for some 7 < fc — 1 and B not divisible by p. This gives 



(3.3) if„ - H, 



n—p*- 



p a BP 



and it follows that the denominator of H n — H n _ p k is divisible exactly by p a . 

We now determine the exponent of the highest power of p dividing (^1), which 
according to Rummer's theorem |12) is the number of borrows involved in sub- 
tracting p k from n in base p. Let n = n; • • • n a • ■ • nf- ■ ■ ■ no be the standard base-p 
representation of n; then n — io = n; • • • ri^O ■ • • 0. Since n Q 7^ and n 3 = for 
k— 1 < j < a, there are a — fc borrows when subtracting p fe from n in base p. There- 
fore the power of p dividing Q^J is p a ~ k , and the power of p in j (iJ„ — H n _ p k ) 
is =p~ k . □ 

Proof of Theorem. Only the terms 1, 2, . . . , n appear as part of the denominators of 
coefficients of the polynomial n\A n (x). Thus, the common denominator is a divisor 
of lcm(l, 2, . . . , n). The previous Lemma shows that every prime power p k < n 
appears. It follows that the denominator of n\A n (x) is lcm(l, 2, . . . , n). 

Since every prime dividing n\ also divides lcm(l, 2, . . . , n), there is no cancellation 
with the numerator of (j)(H n — Hn—j)^ -1 when we divide by nl, so the 

denominator of A n {x) is n! lcm(l, 2, . . . , n). 



Corollary 3.3. The expression /3„ in (|1.7|) is now given by 

(3.4) p n = . y 2 ' 2 '- ,n) iv 

lcm(l, 2, ...,n—l) 

This gives ([T5]) . 



4. The polynomial A n is logconcave 

A sequence of coefficients {ao, 01, . . . ,a n } is unimodal is there is an index 
such that oq < a± < ■ ■ ■ < a,j* and Oj„ > aj*+i > • • • > a n . A polynomial is 
called unimodal if its sequence of coefficients is unimodal. The polynomial is called 
logconcave if its coefficients satisfy > aj_ia.,+i. An elementary argument shows 
that logconcavity implies unimodality |18j . Unimodal and logconcave sequences 
appear frequently in algebra and combinatorics. The reader will find in [51 [TB] a 
survey of these results. 

The logconcavity of the polynomial B n (x) is elementary; it simply corresponds 
to that of the binomial coefficients. The argument for A n (x) is established next. 

Theorem 4.1. The polynomial A n {x) is logconcave. 
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Proof. The result is equivalent to the inequality 



2 



(4.1) II (H n - H n -iY > I . _ J I . + J (F„ - H n - J+1 ) (H n - H^j-x) 
Introduce the notation 

(4.2) / n (j) =H n -H n _ j = ]T — -5 



n 

i=0 

then (|4.1[) is equivalent to 

(4.3) ( " +1)(, ;." j + 1) / n 2 0-)-/n0-) + l>0. 
This is a quadratic inequality with discriminant 

, , j(4n + 5) - 4(1 + n) 2 nUn + 5) - 4(1 + n) 2 3n + 4 

(4.4) — — < — — = < 0, 

3 3 3 

that gives the logconcavity of A n . □ 

Note 4.2. There are many instances where polynomials appearing in connection 
with the evaluation of integrals are logconcave. For instance, the polynomial 

m 

(4.5) P m (a) = ^di^a 1 



m 

k 



with 

(4.6) dt, m = 2- 2m J2 2 

appears in the formula 
(4.7) 



k=l 



,/2m-2fc\ fm + k\ (k 
\ m — k J \ m J \l 



dx 7T Pm(d) 



(x 4 + 2ax 2 + l) m+1 2 



l 



[2(a + l)] ro+ i 

The reader will find in [T] different proofs of these formulas. The unimodality of 
P m (o) was established in [5] and its logconcavity, in [10. . A direct proof of this 
result also appears in [B]. 

The logconcavity of a sequence {aj : 1 < j < n} can be expressed in terms of the 
operator £ defined by £({a,j}) — {a 2 — a,j^ia,j + i}, where aj — if j < or j > n. 
Thus, a positive sequence a is logconcave if £(a) is positive. A sequence is called 
r-logconcave if £' s ^ (a) is positive for < s < r. It is called infinite logconcave if it 
is r-logconcave for any r G N. 

The next conjecture is based on extensive symbolic calculations. 

Conjecture 4.3. For every n £ N, the polynomial A n (x) is infinite logconcave. 

The question of logconcavity of a sequence {afc} is intimately connected with 
the location of the zeros of its generating polynomial P(x) = ciq + a\x + ■ ■ ■ + 
a n x n . Newton showed that if the zeros of P are real and negative, then {a^.} is 
logconcave. It turns out that the zeros of the polynomial P m (a), generated by 
the sequence {dz, m } in (|4.6[) . do not satisfy this condition. The quest for a proof 
of the logconcavity of {d^ m } lead one of the authors to study the operator £. 
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Experimental data suggested that this sequence was infinite logconcave. A similar 
conjecture on the binomial coefficients was proposed as a testing problem. 

The operator £ does not preserve logconcavity. In order to bypass this difficulty, 
McNamara and Sagan |14j introduce a remarkable method that gives infinite logcon- 
cavity. Given r£Ra sequence {a&} is called r-factor logconcave if a 2 > rak+xdk-i- 
The next lemma appears in [14] . 

Lemma 4.4. Let {at} be a non-negative sequence and let ro = (3 + \/5)/2. Then 
{afe} being ro-factor logconcave implies that £({afc}) is too. So in this case {a^} is 
infinite logconcave. 

This result can be used to check that a sequence a n j is infinite logconcave 
for a specific fixed value of n. It was used in [T3] to show that the sequence of 
binomial coefficients {(^) : < k < n} is infinite logconcave for n < 1450. Their 
approach is remarkably simple: calculate £, l {a n j) for i up to some bound M. If 
the sequences £(a nj ), £ 2 (a nj ), . . . , £ M (a nj ) are non-negative and £, M (a n .j) is ro- 
factor logconcave, then Lemma 14.41 implies that {a ra j}™ =0 is infinite logconcave. 
Employing this technique we have verified the following result. 

Theorem 4.5. The polynomial A n {x) is infinite logconcave for all n < 300. 

Extending the conjecture on the infinite logconcavity of the binomial coefficients, 
Stanley [T7], McNamara-Sagan [2] and Fisk [7] conjectured the next result. This 
was established by Petter Branden in [3] . 

n 

Theorem 4.6. Suppose that the polynomial cikX k has only real and negative 

k=0 

zeros. Then so does the polynomial 

n 

(4.8) ^ ( a l - Ofc-iOfc+i) zk , where a_i = a n+1 = 0. 

fc=0 

In particular, the sequence {at ■ < k < n} is infinite logconcave. 

The theorem established the fact that the binomial coefficients { (?) : < k < n} 
are infinite logconcave for every n € N. On the other hand, examples of sequences 
conjectured to be infinite logconcave have generating polynomials with non-real 
zeros. This holds for P m (a) in (|4.5p as well as the polynomial A n (x) considered 
here. 

5. An extension 

The expression for the iterated integrals described here extends to iterated inte- 
grals obtained by 

(5.1) f 0J (x) - In^l + z) 

fn.j(x) = / f n -l,j{t)dt. 
JO 

Symbolic computation gives 

/i >2 (a:) = 2x - 2(x + l)u + (x + l)u 2 

f 2)2 (x) = ^x(7x + 6)-^(x + l) 2 u+^(x + l) 2 u 2 

/s,2(z) = -^-x(85x 2 + U7x + 66)-^(x + l) 3 u+hx + l) 3 u 2 
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where u = ln(l + x). The structure of f n ,j(x) is provided below. 

Theorem 5.1. For n > and j > 1, there are polynomials A n j(x) and B n ^j(x) 
such that 

j 

(5-2) fn,j{x) — A n j(x 

k=l 

Moreover, B n> k,j{x) = b n> k,j(l + x) n , for some 6 n ,fej € Q. 

Proof. In the first step, it is shown that /„ j has the stated form, for some polyno- 
mials Anj, B Ui k t j. This is elementary: integration by parts shows that, for r, k G N, 

(5.3) J s r ln k sds = P(lns)s r+1 , 

where P is a polynomial of degree k. Induction on n and (|5.ip . give the result. 

The next is to prove that B n ^j(x) is a multiple of (1 + x) n . In order to achieve 
this, replace (|5.2I) in (|5.1[) to obtain 

(5.4) (1 + x)A' nJ {x) + B niltj (x) = (l+x)A n ^(x) 

{l + x)B' nkj (x) + {k + l)B n . k+l ,j(x) = (l + x)B n _ l>kd (x) farl<*<j-l, 

B 'nj,j(x) = B n ^ hjJ (x). 

Fix n and j and use induction on k to show B nt k,j(— 1) = 0. The first equation in 
(I5.4|) gives 1) = and the second one gives B nik ,j(—1) = 0, using induction 

on k. The fact that B n ^,j{x) is a power of (x -I- 1)" now follows directly from the 
recurrence (|5.4|) . Indeed, 

1 + x 1 + x 

( 5 - 5 ) S„, fc+1J (x) = — ^Bn-i,fe,i(a;) ~ YTfc- 8 "^^^ 

gives the result for 1 < k < j — 1. The case fc = j is settled with B' n ^{x) — 
B n-ij,j{x), where the vanishing of B n jj(x) at x = —1 adjusts the constant of 
integration. Therefore 

(5.6) B n , kd (x) = b n , k , j {l + x) n 

for some b n ^.j to be determined. 

For n > 1 and j > 1, the coefficients 6 n ,fe,j satisfy the equations 

(5.7) ftojj = 1 

&o,k,j = for < k < j - 1 

k + 1 1 

&n,fcj = b n,k+l,j + -K-l,kj for < fc < j - 1. 

The third recurrence in (15.71) can be solved directly to obtain b n ,j,j — 1/nl. □ 
To write b n>k) j more explicitly, let 

(5.8) H n , m = y 

, , . t-r 1 . . aia 2 ■■■a m 

for n > and m > 0. In other words, the sum is over all nondecreasing m-tuples 
with entries from {1,2,..., n}. H n>m generalizes the harmonic number H n = H n ,x- 
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For m — the product is empty, so i?n.o = 1 for n > 0. For n = the sum is 
empty unless m = 0, so Ho. m — for m > 1. 

For n > 1 and m > 1, we have H n m — —H n ,m-\ + H n _i m by breaking up the 
m-tuples in H n m according to whether a m = n or not. It follows that the proposed 
expression for b n ,k,j satisfies 

k + 1 1 

(5.9) b n h ,j = b n k+i,j H — b n -i k,j- 

n n 

It follows from here that 

(5.10) b n k j = — — H n .j-k 

n\k\ 



satisfies the equations in (|5.7I) . 

Theorem 5.2. The polynomial B n> k,j in Theorem 15. II is given by 
(5.11) B ntktj (x) = tlZ_Jl H n j_ k (l + x) n . 

Note 5.3. The recursion H n m = \ii n ^ m -\ + Hn-i,mi with n fixed, can be solved 
explicitly. The initial term is H± m = 1. It may be checked directly that 



(5-12) H n,m = J2(-l) k+X ( n k )k 



— in 



k=l 

For n — 2, H^^-k is a geometric series; therefore (|5.10p gives 

{-iy- k j\(y- k+1 - 1) 



(5.13) b 2 . kj 



The last step in the determination of f n j is the construction of the polynomial 
A n ,j{x) in (|5~2l) . 

The equation for A n j(x) comes from (|5.4I) : 

(5.14) A' n<j (x) = An-x t j(x) + ( -*] Jj! ff n ,j-i(l + x) n -\ 

with the boundary conditions Aq^(x) — and A n j(0) = 0. 
The polynomial A n ^ (x) is now written as 

n 

(5.15) A n j(x) = S ^ c nr jx r 

for some coefficients c„ irj -. The fact that A n j is of degree n comes directly from 
(15.141) . This recurrence also yields 

(5.16) c n ,\,] = a n ,j 

1 1/n-l 

(5.17) Cn,r,j — l,r— X,j ' 



r V r — 1 

where 
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The claim is that the solution of this recurrence yields the following expression 
for the polynomial A n j(x). 

Theorem 5.4. The polynomial A n j(x) in Theorem 15. II is given by 



r— 1 ^ ' 

Proof. Direct replacement of 



V-l 



E 



Hfi— k,j— l 



fn\ ^ Hn-k,j-i 



(5.20) : E 



n! V r 7 ~* n — k 

v 7 fe=0 

in ([536]) and (|57T7|) . □ 

Note 5.5. The case j = 1 of (|5.19[) reduces to the expression for given 
in Theorem I2TJ 



The final statement deals with an arithmetical property of the coefficients of 
B n j(x). This was first observed symbolically and it was one of the motivations for 
the work discussed here. 

Corollary 5.6. Assume p is a prime that divides the denominator of a coefficient 
of B nt kj{x) in Theorem 15.21 Then p < max(n, k). 
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